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In this work we use geometrical optics and the caustic-touching theorem to study, in an exact way, 
the change in the topology of the image of an object obtained by reflections on an arbitrary smooth 
surface. Since the procedure that we use to compute the image is exactly the same as that used to 
simulate the ideal patterns, referred to as ronchigrams, in the Ronchi test used to test mirrors, we 
remark that the closed loop fringes commonly observed in the ronchigrams when the grating, referred 
to as a Ronchi ruling, is located at the caustic place are due to a disruption of fringes, or, more 
correctly, as disruption of shadows corresponding to the ruling bands. To illustrate our results, we 
assume that the reflecting surface is a spherical mirror and we consider two kinds of objects: circles 
and line segments. 

1. Introduction 

By using the paraxial approximation of geometrical optics, Berry in an extraordi- 
nary work [Ij, among other things, studied the image of an arbitrary one-dimensional 
object obtained by reflection on an arbitrary smooth surface (of rippled water). He 
found that under certain conditions the object and its image do not have the same 
topology. To explain this beautiful phenomenon he introduced the so-called caustic 
touching theorem, which states that changes of image topology occur when the object 
touches the caustic associated with the family of imaginary light rays emitted by the 
observing eye. 

It is worthwhile describing, with little convenient changes, the procedure followed 
by Berry to obtain the image of an arbitrary object curve by reflection on an arbitrary 
smooth surface. Without loss of the generality, we assume that: the object curve is 
lying on a plane perpendicular to the z axis, the two-dimensional reflecting smooth 
surface is locally given by z = f{x,y) and the position of the observing eye by S" = 
(si, 52,53), (see figured]). At first instance (sight), one could think that to obtain the 
image of the one-dimensional object one has to take into account the family of light rays 
associated with each of its points. That is, a family of light rays characterized by three 
parameters (one of these parameters provides the position of an arbitrary point on the 
one-dimensional object and the other two give the direction of the light ray emitted 
from that point). However, the only important light rays that give contribution to the 
image formation are those that reach the observing eye and these, in accordance with 
the reciprocity principle [2], can be regarded as belonging to the single family emitted 
by the eye. In other words, the original problem of image formation is equivalent to a 
new problem where the observing eye is replaced by an imaginary point light source. 



From this new point of view, a point on the reflecting surface belongs to the image 
associated with the one-dimensional object if it can be associated, via a reflected 
light ray emitted by the imaginary point source, with a point of the one-dimensional 
object, as it is graphically described in flgure[2j Since in general the curvature of the 
reflecting surface is not equal to zero, then the reflected light rays emitted by the 
imaginary point light source will focus at a region in the space. This region is the 
caustic associated with the imaginary reflected light rays. From a mathematical point 
of view, the evolution of the reflected light rays is described by a map between two 
subsets of or equivalently by a one-parameter family of maps between two subsets 
of (see section 2). The caustic associated with the reflected light rays is determined 
by looking for the points, on the reflecting surface, where these maps are not locally one 
to one. Therefore, if the one-dimensional object is placed outside the caustic region 
there will be a one to one correspondence between its points and the points of its 
associated image. That is, the observing eye will see only one image. In this case one 
says that the object and its image have the same topology. Remember that, roughly 
speaking, two curves are said to be topologically equivalent if one can be transformed 
into the same shape as the other without connecting or disconnecting any points. 
However, if the object is located at the caustic region the maps are not one to one 
and its position with respect to the caustic will be crucial to describe correctly its 
image. The caustic-touching theorem establishes that new image loops appear when 
the object curve touches the caustic associated with the reflected light rays emitted 
by the imaginary point light source. Such disruption may be elliptic, loop born from 
an isolated point, or hyperbolic, loop pinched off from an already existing one. In 
these cases, the object and its image do not have the same topology, and the observer 
may see several images corresponding to a single object. As remarked by Berry it 
is what happens in nature's optics, whose elements may be the reflecting surface of 
rippled water. The flrst aim of the present work is to obtain, within the geometrical 
optics approximation, an exact set of equations to study the change of topology of a 
one-dimensional object obtained by reflection on an arbitrary smooth surface. 

Prom the research developed by several authors[ll O IU [5l El [71 [8] it is quite clear 
that the multiple image formation process in the natural optical systems, whose com- 
ponents may be the reflecting surface of rippled water, refractive-index gradients in 
the atmosphere, or the gravitational fleld associated with a matter distribution, can 
be explained within the geometrical optics limit by using the caustic-touching theo- 
rem. On the other hand, the conventional optical systems such as microscopes and 
telescopes, whose basic components are mirrors and lenses, are designed to produce a 
single image, of course, as perfect as it is possible. So, at flrst glance, one could infer 
that non-multiple images, associated with a single object, can be seen by using these 
optical devices. However, in general, it is not true. The second and main contribution 
of this work is to remark that even with a perfect mirror it is possible, under certain 
conditions determined by the caustic-touching theorem, to observe image disruption. 
In particular, we remark that the caustic-touching theorem allow us to describe the 
pattern, referred to as the ronchigram in the well known Ronchi test, when the grating 
is placed at the caustic region. 

The Ronchi test [9] developed by Ronchi in the 1920's is one of the most simple and 
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Figura 1: Schematic drawing of the object curve, which we assume is lying on a plane perpendicular 
to the z axis, the two-dimensional reflecting smooth surface locally given by z = /(a;, y] and the 
position of the observing eye given by S = (si, S2, S3). 



powerful methods to extract information about the aberrations of an optical system. 
For these reasons it has been the subject of innumerable publications, from both the 
physical and the geometrical points of view (see Ref. jlOj and the references cited the- 
rein) . From the physical point of view the fringes are interpreted as due to interference 
between several wavefronts produced by the ruling acting as a diffraction grating, while 
from the geometrical point of view the fringes are interpreted as shadows of the ruling 
slits. When the frequency of the ruling is not very high, the two points of view predict 
the same result. Since by analyzing the fringes of the real and ideal ronchigrams, one 
can determine the type and, in principle, the magnitude of the aberrations present 
at the exit pupil of the system, then it is fundamental to know the properties of the 
ronchigrams. 

In this work we consider the geometrical point of view of the Ronchi test. Under 
this assumption, the essential features of the Ronchi test for a concave mirror when the 
point source is located on the optical axis may be described by reference to figure[3l 
The light rays emitted by the point light source are reflected by the mirror under test 
and they focus to a region in the space. As previously mentioned, this region is the 
caustic associated with the reflected light rays. The grating, which is referred to as 
the Ronchi ruling, is located at different positions on the optical axis. The pattern 
observed through the grating on the surface of the mirror is referred to as the real 
ronchigram. From the geometrical point of view, the fringes of the real ronchigram 
are interpreted as shadows of the ruling bands. By comparing these real fringes with 
the ideal ones obtained by simulation one can deduce the defects of the mirror under 
test|llj. Therefore, one of the main corner stones of the Ronchi method is the ideal 
ronchigram which is obtained by simulation. Thus, Sherwood [12] has calculated the 
Ronchi pattern of a paraboloidal mirror when tested near the center of curvature; 
Malacara developed an algorithm to predict the geometrical ronchigram of any 
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Fi gura 2: The new geometrical arrange to compute the image of a one-dimensional object un- 
der reflection on the arbitrary smooth surface z = f{x,y). The observing eye has been replaced by 
an imaginary point light source. In this diagram we show an imaginary emitted light ray and the 
corresponding reflected light ray which, we assume arrives at a point of the one-dimensional object. 
Therefore, the point of the surface where the imaginary light ray is reflected belongs to the image of 
the one- dimensional object. 



spherical or aspherical mirror when tested at any point along the optical axis. These 
results have been generalized by Cordero, Cornejo and Cardona [13j to an arbitrary 
mirror when the point source is located at any position, in particular their equations 
can be used to simulate ronchigrams for the cases of centered and off-axis conic sections 
with the point light source at any location. Cordero, Diaz, and Cabrera [Hj have 
presented a simple algorithm that allows the simulation of ronchigrams for any optical 
system in which it is possible to make an exact ray tracing. In accordance with the 
results obtained by these authors it is clear that the form and structure of the ideal 
ronchigram associated with a given perfect mirror depend on the location of both the 
point source and the Ronchi ruling. Furthermore, it is well known that when the Ronchi 
ruling is placed at the caustic region it is common to see closed loop fringes in the 
ronchigram. However, to our knowledge, an explanation of the closed loop fringes in the 
ronchigram using geometrical optics has not been presented. Actually, in practice for 
an easy interpretation of the ronchigram it is common to avoid any closed loop fringe, 
which is possible only when the grating is located outside the caustic region associated 
with the reflected light rays [lOj. In this work we remark that the explanation is given 
by the caustic-touching theorem, and therefore, the relative position of the grating with 
respect to the caustic associated with the reflected light rays emitted by the point light 
source is crucial to describe correctly the structure of the associated ronchigram. Our 
justification for these assertions is as follows. From figures [2] and [3l it is clearly obvious 
that the procedure followed by Berry to compute the image of an arbitrary object 
under reflection on an arbitrary smooth surface when the observer is located at an 
arbitrary position in the space is exactly that used by several authors to simulate 
the ronchigram associated with an arbitrary smooth reflector by using a point light 
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source located at different positions in the space. Remember tliat in Berry's procedure 
the observer is replaced by an imaginary point light source, if in addition, we replace 
the one-dimensional object by a grating, in particular a Ronchi ruling, then we get 
the desired configuration to simulate the corresponding ronchigram. Even though the 
simulation of the ronchigrams has been implemented for many years, to our knowledge, 
nobody has remarked on the fundamental role played by the caustic associated with 
the reflected light rays to describe the structure of the ronchigrams when the Ronchi 
ruling is located at the caustic region. Actually, as we will show in section 2, the 
equations that we obtain to compute the image associated with an arbitrary object 
obtained by reflection on an arbitrary smooth surface, are exactly those obtained by 
Cordero, Cornejo and Cardona [13j in the context of the Ronchi and Hartmann tests. 
However, these authors did not study the properties of the ronchigrams when the 
Ronchi ruling is located at the caustic place. 




Figura 3: Schematic drawing of the Ronchi test arrangement. In this diagram we have the surface 
under test, locally given hy z — f{x,y), a real point light source located on the optical system and a 
Ronchi ruling. The pattern observed through the grating on the surface of the mirror is referred to as 
the real ronchigram. 

An important observation is that in the Ronchi test the point light source is not an 
imaginary one, as in Berry's procedure, but it is a real one and the pattern, the image 
(in Berry's problem), is observed through the grating on the surface of the mirror. 
For practical purposes it is more convenient to plot the simulated ronchigram in a 
plane close to the reflecting surface. Therefore, in this work we plot the corresponding 
images (in Berry's problem) or patterns (in Ronchi's test) in the plane z = 0. 

The aim of the present work is twofold: first, within the geometrical optics limit, 
we obtain an exact set of equations to study the changes of image topology under 
reflection on an arbitrary smooth surface and, second, we apply our results to two 
particular cases to show that even with a perfect spherical mirror one can observe 
the change of topology, in particular we remark that the structure of the ronchigrams 
when the Ronchi ruling is located at the caustic region can be described by using the 
caustic-touching theorem. 



5 



The organization of the present work is as follows: In section 2, we use geometrical 
optics to obtain an exact set of equations to compute the image of a one-dimensional 
object produced by reflection on an arbitrary smooth surface. We find that this set of 
equations is exactly that previously obtained in the context of the Ronchi and Hart- 
mann tests in[T3]. In section 3, we compute the caustic surface associated with the 
light rays reflected by an arbitrary smooth surface when a point light source is located 
at an arbitrary position in the space. In section 4, we collect the two sets of equa- 
tions to study in an exact way the changes of image topology. Finally, in section 5, we 
present some examples to illustrate the change of image topology when the reflecting 
surface is a spherical mirror and the objects are circles or line segments. In particular, 
we describe the structure of the ronchigram associated with a spherical mirror when 
the point light source is located on the optical axis and the Ronchi ruling is placed at 
the caustic region. 

2. Computation of the image of a one-dimensional object 
obtained by reflection on an arbitrary smooth surface 

In this section, following Berry's procedure, we obtain the exact set of equations to 
compute the image of an arbitrary one-dimensional object obtained under reflection on 
an arbitrary smooth surface. Remember that the original problem is changed by a new 
one, where the observer is replaced by an imaginary point light source as previously 
explained and shown in figure[2j Since this new problem is exactly that of simulating 
ronchigrams in the Ronchi test and there the point light source is a real one, from 
now on, we will use the term point light source to refer to both the real and imaginary 
point light sources, understanding that in the Ronchi test it is real and in the Berry's 
procedure is an imaginary one. Thus assuming that in TZ^ we have a smooth arbitrary 
surface and a point light source with arbitrary position we first obtain a parametric 
representation of the maps that describe the evolution of the reflected light rays by 
the arbitrary smooth surface, and then we obtain the set of equations that allow us to 
compute the image of a one-dimensional object obtained by reflection. 

If the point light source is located at S = {81,82,83) and the smooth arbitrary 
surface is locally given hy z = f{x,y), see figure HI then the light ray reflected at the 
point r = {x,y, f{x,y)) is described by 



is the unit normal vector to the reflecting surface z = f{x,y), and / gives the 
direction of the diverging ray from the point light source. Prom figured] we have 



T = r + lR 




where / is the distance along the reflected light ray, R is given by 



R = i -2{i ■N)N, 



(2) 





(3) 



X-8^Y + {y-82Y + {f-83Y 
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Figura 4: Schematic drawing of the optical system and the vectors used to compute the image 
of a one-dimensional object obtained by reflection on an arbitrary smooth surface locally given by 
z = f{x,y). In this diagram S — (si, 52,53) denotes the position of the point light source, / the 
direction of an emitted light ray, f = {x, y, f(x, y)) the point on the smooth surface where the emitted 
light ray is reflected in the direction R and N is the normal vector to the smooth surface at the point 
of reflection. 



In order to obtain a vector field perpendicular to the reflecting surface, we define 
the function G{x,y,z) = z — f{x,y). Observe that one level surface of this function 
is the reflecting surface z = f{x,y). Therefore, a vector field perpendicular to the 
reflecting surface is given by 

N={-f,,-fy,l), (4) 

where fx = (df /dx)\G=o and fy = {df /dy)\G=Q- Finally, the unit normal vector field 
to the reflecting surface is given by 

N = -thzML. (5) 



By using Eqs. (l)-(5), a direct computation shows that if 5 = (81,82,83) is the 
position of the point light source then a light ray that is emitted in the direction / 
and reflected by the arbitrary smooth curved reflector at the point r = {x,y, f{x,y)) 
is described by 

Ti = x + (6) 
a 

T. = , + ^, (7) 

a 

n = f{x,y) + ^, (8) 
a 
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where 

hi = {x - si){l - + f^) - 2Mfy{y - S2) + S3 - f], 

h2 = {y-S2){l + f!-f^)-2fy[Ux-si) + S3-f], 

hs = {f-S3){-l + f^ + f^) + 2[Ux-si) + fyiy-S2)], 

a = (1 + f! + fyWisi - xY + {S2 - yf + {S3 - (9) 

and /, which gives the position of an arbitrary point along the reflected light ray, is 
such that f{x,y) < I < oo. Furthermore, we assume that Xmin < x < Xmax and 
Vmin < y < Vmax- The valucs of x^in, Xmax, Vmin and ymax are determined by the 
dimensions of the reflecting surface, which we assume known. 

Since we are interested in the intersection of the reflected light rays with an arbi- 
trary plane z = constant, then it is convenient to take T3 = zq, and thus from Eq. ([8]) 
we have that 

, f T3- fix,y) \ f Zq- fix,y) \ 

so that Eqs. ©-(HI) can be rewritten in the following form 

Ti{x,y,zo) = x + [zo- f{x,y)] | ^' '^2, .§3) j ^ 

rj. ( N , r ft M ( h2{x,y,Sl,S2,S3) \ 

T2{x,y,zo) = y+[zo- f{x,y)\[—, r , 

\h3ix,y,Sl,S2,S3)J 

T3{x,y,zo) = zo, (11) 

where f{x,y) < zq < 00. 

Before continuing it is important to explain the geometrical meaning of the pa- 
rametric map given by Eqs. (jlip . For this end, we assume that the position of the 
point light source is fixed. Then for fixed values of x and y; that is, for a point on the 
reflecting surface, z = f{x,y), Eqs. (jlip describe a line segment, which starts at the 
point {x, y, f{x, y)) and goes to infinity as zq does. Therefore, as x and y take all their 
allowed values, this map describes a family of line segments that start at the points of 
the reflecting surface and go to infinity. This family of line segments are the reflected 
light rays. On the other hand, from a mathematical point of view, Eqs. (jlip describe 
the parametric form of a map between two subsets of Ti?, where (x, y, zq) are local 
coordinates of the domain space and (Ti, T21 T3) are local coordinates of the target 
space. Equivalently, Eqs. (jlip . can be seen as a one-parameter family of maps between 
subsets of Ti?, where each map is characterized by a particular value of zq. This family 
is explicitly given by 

Ti{x,y,ZQ) = x + [zo- f{x,y)] | ^' ^i' ^2, ^3) j ^ 

\h3{x,y,Si,S2,S3)J 

rr ( N , r f( ^1 ( h2{x,y,si,S2,S3)\ 

T2(x,y,zo = y+[zQ-f{x,y)\[— r . 12 

\h3[x,y,Si,S2,S3) J 

Each member of the family, characterized by a specific value of zq, maps points of 
the reflecting surface to points on the plane z = zq. It is important to remark that in 
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Eqs. (|12p . zq is not a coordinate as it is in Eqs. (jlip . but it is a parameter characterizing 
a particular member of the family of maps, and is such that f{x,y) < zq < oo. 

The sets of Eqs. (jlip and (jl2p are those maps we were referring to in the introduc- 
tion, which describe the evolution of the reflected light rays by the arbitrary smooth 
surface after being emitted by the point light source. We remark that because we have 
used the reflection law only one time to obtain Eqs. (jlip. and, therefore, Eqs. (jl2p . 
these equations describe the evolution of the light rays that have experienced only one 
reflection before leaving the smooth surface. In this work we assume that the parame- 
ters that characterize the optical system under study are such that this condition is 
fulfilled by the reflected light rays. 

Assuming that in the plane z = zq a coordinate system (T^, Ty) is introduced, with 
origin at (0, 0, zq) such that and Ty are parallel to the axes x and y respectively (see 
figure[5]), a one-dimensional object lying on this plane can be described in a parametric 
way by 

= r(a), 

Ty = S(a), (13) 

where o" is a parameter that labels the points on the object. If we eliminate the para- 
meter a we have that the object could be described by 

T. = MTy). (14) 

Therefore, the image of this object that an observer, with position S = {31,82,83), 
may see on the surface of reflection is given by all the points of the form (x, y, f{x, y)) 
such that X and y are solutions to 



T^[a) = x+[zo- f{x,y)] 



hi{x,y, 81,82,83] 

h3ix,y,8i,82,83] 



^1 \ , r~ ft ^^ ( h2{x,y, 81,82, 83)\ 

Tyicy) = y + [zo-f{x,y)][— . 15 

\h3[x,y,8i,82,83) J 

In other words, a point on the reflecting surface belongs to the image of the object 
(jl3p if it can be connected, via a reflected light ray, with a point of that object (see 
figure Ej) . In the explicit examples that we present later on, we will not plot the image 
on the reflecting surface, but on the plane z = 0. This is because when one takes a 
photograph of an object its image is normally printed on a plane. Therefore, in our 
examples, the image will be given by all the points {x, y, 0) such that x and y satisfy 
Eqs.([l5]). 

Equations ()15p together with the caustic associated with the reflected light rays 
will allow us to study the change in the topology of the image of an arbitrary one- 
dimensional object obtained under reflection on an arbitrary smooth surface. 



3. Computation of the caustic 

Since to study the change of image topology of an arbitrary one-dimensional object 
obtained by reflection it is required to compute the caustic associated with the reflected 
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Figura 5: Schematic drawing of the two sets of coordinate systems used to compute the image 
of a one-dimensional object obtained by reflection on an arbitrary smooth surface locally given by 
z = f{x,y). We also have included an emitted light ray such that its associated reflected light ray 
connects a point of the smooth surface with a point of the one-dimensional object. The point, on the 
smooth surface, where the emitted light ray is reflected belongs to the image of the one- dimensional 
object. 



light rays described by Eqs. ([IT]) or equivalently by Eqs. (fT2|) . in this section, following 
our previous work [HIITB], we review the computation of the caustic by using Eqs. (fTTj) 
and furthermore we point out its relationship to the caustic associated with Eqs. (fT2]) . 
To this end, we introduce the following: 

Definition: Let /i : — > be a differentiable map, with A4 and M differentiable 
manifolds. The set of points in Ai where h is not locally one to one are referred to as 
its critical set, and the image of the critical set is referred to as the caustic set of h 
[n\ [T8l I19j . If M and M are differentiable submanifolds of TZ"" with local coordinates 
(xi) and {yj) respectively, then locally h is given by 

yi = hi{xj), where i,j = l,...,n. (16) 

Therefore, in this case the critical set is obtained from the condition 

..de.(|i)=0^ (17) 

In general this condition can be written in the following way 

F{xu...,xn) = 0, (18) 

if this equation can be solved, for example, for Xn, then locally the critical set of h is 
given by 

Xn = g{xi,...,Xn-i), (19) 

which, in the general case, is a set of surfaces of dimension n — 1 in the domain space 
which we are assuming has dimension n. Therefore, the caustic set, which is the image 
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of the critical set, is obtained substituting Eq. (fT9|) into Eqs. p^ . That is, the caustic 
set associated with the map ()16p is locally given by 



yi = hi{xi,...,Xn-l,g{xi, ...,Xn-l)), 
2/2 = h2{xi,...,Xn-l,g{xi,...,Xn-l)), 



Vn 



hn{xi, ...,Xn-l,g{xi, ...,Xn-l))- 



(20) 



In accordance with the above definition, the critical set of the map given by 
Eqs. (jlip . that is, the set of points in the domain space with coordinates (x, y, zq) 
such that the map is not locally one to one, is obtained from the following condition 



J{x,y,zo) = det 



diTi,T2,T3) 



\ dx J \ dy J \ dy J \ dx 



0. (21) 



d{x,y,zo) 

By using Eqs. (jlip a direct computation shows that this condition is equivalent to 

2 



J{x,y,zo) = H2{x,y) 



zq- f 



+ Hi{x,y) 



+ Ho{x,y) =0, 



(22) 



where 



H2{x,y) 

Hi{x,y) 
Ho{x,y) 



h- 
h- 



dh\ dh 
dx j \ dy 



f dr\ 

\dx J \dy 



dx J \dy 



dh 
dy 

df 



( dh\ f df 
+ 7^ X — 



dx 



\dy 



(23) 



with 



r 
h 



ix,yJix,y)), 

(/ll,/l2,/l3)- 



(24) 



From Eq. (j22p we find that the critical set associated with the map given by Eqs. ([TT 
is given by 

^ -Hi ± ^Hf-4H2Ho\ 



Zq = zo±{x,y) = f + h3 



2Ho 



(25) 



Therefore, the caustic set associated with the map given by Eqs. (jlip . which by defini- 
tion is obtained by substituting Eq. (j25|) into Eqs. pT]) . can be written in the following 
form 



^^^=^+ r- 



(26) 
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It is important to remark that this equation is equivalent to that obtained by Shealy 
and Burkhard |20il21j by using a different procedure. From this last equation it is clear 
that the caustic set or simply the caustic associated with the reflected light rays des- 
cribed by Eqs. (jlip . in general, is composed by two branches, which for very particular 
forms of the reflecting surface and particular positions of the point light source reduce 
to a single point. For example, this happens when the reflecting surface is part of a 
perfect spherical mirror and the point light source is located at its center of curvatu- 
re. In the general case the two branches of the caustic are two-dimensional surfaces 
which when they are stable under small deformations of the reflecting surface and the 
position of the point light source, they locally have singularities of well known types: 
the swallowtail, the pyramid or elliptic umbilic and the purse or hyperbolic umbilic 
[T71 [T8t ri9] . The fact that the caustic might not be stable under small deformations 
of the reflecting surface and the position of the point light source is related to the 
symmetries of the system formed by the reflecting surface and the point light source. 

Since we are assuming that Xmin ^ 2; < x^ax and ymin ^ y ^ Umax 7 then in our 
case the caustics will be located within a finite region of the space. In particular there 
exist real numbers Zcmin± and Zcmax± such that Zcmin± < T^c± 

Now we compute the critical and caustic sets associated with the one-parameter 
family of maps given by Eqs. ()12p and we point out their relationship to those asso- 
ciated with the map (fTT|) . To this end, we select the map given by Eqs. (fT2|) when 
= zq = constant. That is, 

V/i3(a;,y,si,S2,S3)y 

rr I \ , \~ f( ( h2{x,y,si,S2,s^)\ 

T2{x,y) = y + [zo-f{x,y)][j—, r . 27 

\hs{x,y,si,S2,S3)J 

where zq is a real number such that f{x, y) < zq < 00. Remember that this map sends 
points from the reflecting surface to points on the plane zq = zq = constant. The 
critical set associated with this map is obtained from the following condition 

or equivalently from 

J>, y) = H2{x, y) (^^) ' + Hi{x, y) {^^) + ^o(x, y) = 0. (29) 

In this equation zq is a number, it is not a variable as zq is in Eq. (j22p . Observe that 
the relationship between J{x,y) and J{x,y,ZQ) is given by J{x,y) = J{x,y,ZQ). This 
means that only when 

Zcmini: ^ Zq ^ Zcmaxzsii ('^^) 

the map given by Eqs. (j27p will be not locally one to one. Therefore, if zq satisfies the 
above condition, from Eq. ()29p we have that locally the critical set associated with the 
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map (|27p may be written in the form 



y = -^{x,zq). 



(31) 



In this last equation we have included zq; which is a number, only to show that it is 
the critical set associated with the map given by Eqs. (|27|) . The caustic set associated 
with this particular map, which is obtained by substituting Eq. ([3T|) into Eq. ([27|) . is 
given by 



in the general case, it is a curve or family of curves in the plane zq = zq = constant. 
For very particular cases it reduces to a point. When the caustic is stable under small 
deformations of the optical system formed by the reflecting surface and the point light 
source, then in accordance with the Whitney's singularity theory |17 t ll8 t [T9]. it locally 
has singularities of fold or cusp type. 

Prom the computations presented in this section it is clear that the caustic asso- 
ciated with the map ([27]) given by Eq. ([32|) is equal to the intersection of the caustic 
associated with the map (fTT]l . given by Eq. ([26ll . with the plane zq = zq = constant. 

Before closing this section, we explain the geometrical meaning of the caustic as- 
sociated with the evolution of the reflected light rays by the smooth arbitrary sur- 
face when the point source is located at any position of the space. To this end, 
consider the pencil of light rays reflected by the differential surface dxdy of the re- 
flector. Prom Eqs. (j27p we have that when zq = f the cross-sectional area of this 
pencil of rays is exactly dxdy, and as the light rays evolve, this area is given by 
dTidT2 =1 J{x,y) \ dxdy =\ J{x,y,ZQ) \ dxdy. If (x, y) belongs to the critical set 
of the map given by Eqs. (p7j) then J{x,y) = J(x,y,io) = and therefore in that 
case the cross-sectional area of the pencil of rays collapses to zero. This result shows 
that the caustic is defined by the focusing region associated with the reflected light 
rays and therefore also can be defined as the singularities of the flux density [20, ^21] . 
Pinally, we remark that in the literature one can find explicit expressions and plots 
of the caustic associated with particular refiectors and certain positions of the point 
light source [22| [231 [21]. 

4. Equations to study the change of image topology 

In this section we summarize the two sets of equations that allow us to study, in an 
exact way, the changes in the topology of the image, of an arbitrary object, obtained 
by reflections on an arbitrary smooth surface: 

Jft The image of an arbitrary one-dimensional object, locally given by Eqs. (113p and 
located on the plane z = zq = constant, obtained by reflection on the smooth reflecting 




(32) 
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surface locally given by z = /(x, y), is computed by solving the Eqs. p^ : that is 



for X and y. The solution must be such that Xmin < x < Xmax and ymin < y < Vmax- 

XX The caustic associated with the reflected light rays described by the map ([271) 
given by the intersection of the caustic (f26l) . that is 



and the plane z = zq = constant. Or equivalently given by Eq. (|32p. 

Remember that the map given by Eqs. (j27|) maps points on the reflecting surface to 
points on the plane z = zq = constant. If the plane z = zq = constant, which contains 
the object, is outside the caustic region; that is, there is not intersection between this 
plane and the caustic given by Eq. ()34p . then the map ()27p is locally one to one. This 
means that for each point {Tx{a), Ty{a), zq) on the object there is a unique solution (x, 
y) to Eqs. ()33l) and therefore, under this condition, the observer will register a single 
image. In other words, the object and its associated image, in this case, have the same 
topology. Now we assume that the plane z = zq = constant is located at the caustic 
region given by Eq. ()34p . In this second case, the map given by Eqs. p7p is not locally 
one to one at those points that belong to the intersection of the caustic and the plane 
z = Zq = constant. Therefore, if the object lying on this plane is outside the caustic 
then the observer will register a single image, but if the object reaches the caustic in 
such a way that they become tangent to each other; that is, there is a touch between 
them, then the object and its associated image do not have the same topology and 
therefore the observer can see multiple images corresponding to a single object (see 
the examples in the next section). 

We close this section with the following observations: 

X) The set of Eqs. ()33p have been reported in the literature by Cordero, Cornejo and 
Cardona [T3] in the context of the Ronchi and Hartmann tests. These authors showed 
that Eqs. ()33p allow to describe the main features of both tests. In particular, they 
found that the ronchigram associated with an arbitrary reflecting surface, when the 
point source is located at an arbitrary position and the Ronchi ruling is at the plane 
z = Zq = constant with its rulings parallel to the x axis, is given by the level curves 



XX) The caustic, Eq. ([M]) . has also been reported in [161 ISOl IH] and, in particular, 
it was used to compute the circle of least confusion associated with a rotationally 




(33) 




(34) 



oiTy{a). 
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symmetric mirror when the point hght source is located on the optical axis |15l \W[ [211 

mi- 

4^44) Though the two sets of equations ([33j) and ([M]) have been reported, to our 
knowledge, they have not been used to explain the structure of the ronchrigrams when 
the Ronchi ruling is located at the caustic region. The main contribution of this work 
is to realize that the procedure implemented by Berry is equivalent to that used in 
the Ronchi test and therefore that the closed loops fringes observed in the Ronchi 
pattern when the Ronchi ruling is located at the caustic region can be explained by 
using the caustic touching theorem. (The present work can be considered as a logical 
continuation of research reported in Refs. [131 (13 IS].) 



5. Examples: Spherical mirror with the point light source 
on the optical axis 

In order to illustrate everything we have presented in the previous sections we 
assume that the reflecting smooth surface is a part of a spherical mirror with radius r 
and diameter D given by 



z = f{x, y) = r — \/r2 — — (35) 

where —D/2 < x < D/2 and —D/2 < y < D/2. Equivalently, this equation can be 
written in the following form 



z = f{p) = r - \/r2 - (36) 
where p = -^/x^ + y2 jg g^ch that < p < Pmax = D/2. For this case 



Ti{x,y,zo) = x[l + {zo-r+y /r^-p^ )G{p)], 
T2{x,y,zo) = y[l + {zo - r + r'^ - p^)G{p)], 
T2{x,y,zo) = zo, (37) 



where 



^(^) ^ 2(r-g)(r^ -p^) -rVr2 -p2 ^^^^ 

\Jr'^ — p^\2p^{r — s) + r2(s — r + \Jr'^ — p^)] 

A direct computation shows that the two branches of the critical set associated with 
the map Eq. ([57]) are given by 



r\2r^ + 2s/9 — 2r^(s + y/r^ — p^) - r{2p^ — s^J — p^)] 

ZQ- = . : 

2^3 _ 2rp2 j_ 2sp^ - r2(2s + Vr2 - p^) 



r^(4r^ — 5rs + 2s^) — (r — s)[r^(4r — s) + 2(r - s)p'^]yjr'^ - p2 

■^0+ ~ = , 

r2 [3r2 — 4rs + 2^2 + 3(s — r)\/r'^ — p^\ 



(39) 
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and the corresponding branches of the caustic set by 



and 



Tic-{x,y) = 0, 
T2c-ix,y) = 0, 
T3c-{x,y) = Z0-, 



Tic+{x,y) 



Tic+ix,y) 



2x{r — s) p 



2„2 



2y{r - sfp^ 



^2J3j,2 _ 4j,g _j_ 2^2 _|_ 3|'g _ r) Y^r^ — p^] 
T-ic+{x,y) = ZQ+. 



(40) 



(41) 





Figura 6: Intersection of the caustic given by Eqs. 
case r = 2415 mm, D = 1470 mm = 2pmax, z = zo 
at (0, 0, 1350 mm). 



20} and (|41|l with the plane y — 0, for the special 
= 11000 mm and the point light source is located 



Equations (j40p describe a hne segment on the optical axis, that is on the z axis. 
This part of the caustic corresponds to the fact that the rays that leave the points in 
a circle on the mirror with its center on the optical axis converge at a point. All the 
convergent points of all the rays that leave all the possible circles on the mirror with 
centers on the optical axis form this part of the caustic (the line segment). By contrast 
Eqs. ()4ip . in general, describe a surface of revolution with a degenerated singularity 
of cusp type. If we intersect this surface with a plane containing the z axis then we 
obtain a curve with a singularity of cusp type. However, if we intersect it with a plane 
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Figura 7: Intersection of the two branches of the caustic with the plane z = 11000. 



perpendicular to the z axis then we obtain a circle of radius yjT'^^j^ + Tl^^- It is an 
example where one of the branches of the caustic is not a two-dimensional surface, this 
fact is related to the axial symmetry of the optical system under study. Even more, 
as it is clear from Eqs. (I39p -(j4ip for the very particular case s = r, that is, when the 
point light source is located at the center of the spherical mirror, the two branches of 
the caustic reduce to an isolated point. 

By using Eqs. (j39|) . a direct computation shows that when s>r then 

2s — r 

T S 

Tsc+iPmax) < Tsc+ip) < • (42) 

2s — r 

Therefore, there is an interval on the z axis where T^^-ip) = ^3c+(p) (see figureE]). This 
means that, in this case, the intersection of the caustic with a plane z = zq = constant 
such that zq = T^c-{p) = T3c+(/o); is a circle and an isolated point (see figureE]). The 
circle is obtained from the intersection between the plane and that part of the caustic 
that is a surface of revolution, while the point is obtained from the intersection between 
the plane and the line segment corresponding to the other branch of the caustic. 

For this case, the one-parameter family of maps between points on the spherical 
mirror and points of an arbitrary plane z = zq = constant, is explicitly given by: 

Ti{x,y) = x[l + {zQ-r + yjr"^ - p^)G{p)l 

T2{x,y) = y[l + {zo-r + y^r^-p^)G{p)]. (43) 

From the previous discussion we have, if zq is such that zq = T^c-{p) = ^3c+(p) then 
the caustic associated with this map is a circle and a point which coincides with the 
center of the circle. 

If in the plane z = zq = constant, we have a one-dimensional object given by (jlSI) . 
then its image is obtained by solving, for x and y, the following set of equations: 

T,{a) = x[l + (zo-r + Vr2-p2)G'(/,)], 

Ty{a) = y[l + (zo-r + v^;^^)G(p)]. (44) 
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Given Tx{cr), Ty{a), zq ^-^d r, the image on the spherical mirror is given by all the 
points (x, y, r— y^r^ — — y^), where x and y are solutions to Eqs. (j44p . In this work, 
instead of plotting the image on the smooth reflecting surface, we will plot it on the 
z = plane. 




Figura 8: (a) Object space and (b) Image space for k = —7. 
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(a) (b) 



Figura 9: (a) Object space and (b) Image space for k — —Rc. 

For the explicit computations we take: r = 2415 mm, D = 1470 mm = 2pmax, the 
point light source at (0, 0, 1350 mm) and z = zq = 11000 mm as the plane where 
the object is located at. Using these data, in figure (l6|), we present the intersection 
of the two branches of the caustic given by Eqs. (140p and ()4ip with the plane y = 0. 
While in figure ([7|) we present the intersection of these two branches of the caustic 
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with the plane z = 11000, that is, with the plane where the one-dimensional object is 
located at. The intersection is a circle of radius Rc = 4,1643 mm and its center. (In 
the plots presented in the next section corresponding to the object space, the plane 
z = 11000, we only show that part of the caustic that corresponds to the circle, we 
are not including the part corresponding to the isolated point.) 

5.1. First example: Circular Object 

In this example we assume that the object is a circle of radius 1 given by 

(T, - kf + = 1, (45) 

where /c is a real constant. That is, we are interested in obtaining the image of the 
object that is a circle with radius equal to 1 with center at (A;, 0), which is located in 
the plane z = 11000. From Eqs. (I44p and ()45p we obtain that for a fixed value of fc, the 
image of the circle (j45p in the plane z = is given by all the values of x and y such 
that: 

{x[l + (11000 - r + ^/r"^ - p^)G{p)\ - kf 

[1 + (11000 -r + y72r^)G(p)]2 = 1. (46) 




y T ^ g -400 -200 200 400 

(a) '« 



Figura 10: (a) Object space and (b) Image space for k = —1,5. 



We have written a computer program in Mathematica to study the change of 
topology of the image for this case and those cases presented later on. That is, we 
give a particular value to k, and in the object space (T^., Ty, 11000) we plot the object 
(circle of radius 1) and the caustic (remember that we only will plot the circle of 
radius Rc). After that, the set of Eqs. ()46|) are solved for x and y under the conditions 
-Pmax = -735 < X < pmax = 735 and -pmax = "735 < y < Pmax = 735. Finally, 
these values obtained for x and y are plotted in the plane z = 0. In what follows, we 
present some plots where the change of topology is clearly observed. 
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In figure ([8|) we show the object and image for A; = — 7. In this case, we only get 
one image because the map that sends points from the spherical mirror to points of 
the object is one to one; that is, the object is out of the caustic region. In accordance 
with the caustic touching theorem we will observe a change of image topology when 
the object and the caustic touch each other, that is when they become tangent to 
each other. If we assign values to k from —7 to 7, then the first change of topology is 
obtained when k = —Rc — 1 (remember that Rc = 4,1643). For this value of A; a new 
image appears from an isolated point. In figure ([9]) we present the case corresponding 
to k = —Rc, here it is clearly observed the second image. There is another change of 
topology when k = —Rc + 1, for this value of /c a third image appears. In figure (llOp 
we present the case k = —1,5, where the third image is clearly observed. Other change 
of topology is obtained for k = —1 because the caustic, in the plane z = 11000, is 
constituted by the circle of radius Rc and its center. For this particular value of k the 
object is tangent to that part of the caustic that corresponds to the isolated point. 
For k = 0, the images are three concentric circles (see figure [TT]) . Now it is clear that 
there will be other changes of image topology when k = 1, Rc — 1 and Rc + 1. 

Finally, in figure ()12p we show a set of circles with their centers on the axis in the 
object space and in figure ([13]) we present the corresponding images. While in figure 
()14p we show a family of objects corresponding to concentric circles with different 
positions with respect to the caustic, and in figure (jlSp we present their associated 
images. It is important to remark that Murty and Shoemarker [26] have presented the 
theory of a method to test optical systems, which is similar to that of Ronchi, but 
instead of using straight lines for the grating they used concentric circles. They have 
presented characteristic patterns for the usual aberrations of optical systems. In some 
of these patterns it is possible to see the change of image topology. However this fact 
was not pointed out by these authors. 
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Figura 12: Set of circles with their centers on 
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Figura 14: Set of concentric circles with their centers on T, 




Figura 15: Image of the set of concentric circles with their centers on Tx- 
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5.2. Second example: Linear object 




Figura 16: (a) Object space: set of line segments parallel to the Ty axis with Ty £ [—6, 6] and the 
caustic, which is a circle of radius Rc = 4,1643mm and its center, (b) Image space: the corresponding 
images. In the Ronchi test the set of lines (a) is the grating or Ronchi ruling and its image (b) is 
referred to as the associated Ronchigram. 



In this second example, which we beheve is the most relevant due to its connection 
to the famous Ronchi test used to test in particular conic mirrors, we assume that the 
object is a line segment parallel to the Ty axis. Actually in the Ronchi test the grating 
is formed by ruling bands, for our purposes we can assume that the line segments that 
we consider, coincide with the center of the ruling bands of the Ronchi ruling. That 
is, we assume that the object is given by 

= k, (47) 

where A: is a real constant. In other words, we are interested in obtaining the level 
curves associated with the function Ti, given in Eqs. ([l3j) . From Eqs. ([Hj) and we 
have that for k fixed, the image in the plane z = is given by all the points of the 
form (x, y, 0) such that x and y satisfy the following equation 

x[l + (11000 - r + \/r2 - p'^)G[p)] = k. (48) 

In exactly the same way as in the previous case, to study the change in the topology 
of the image, we give different values to k in the interval [—6, 6]. In figurefTBTa) we 
show the caustic and the object for different values of k. In figurefTUlfb) can be observed 
the corresponding images. In the Ronchi test figurefTBTa) corresponds to the Ronchi 
ruling and figurefTUTb) corresponds to what is referred to as the ideal ronchigram 
associated with a spherical mirror. The ronchigram, figurefTBTb) . has been reported 
in the literature, but its closed loops fringes never had been explained by using the 
caustic touching theorem as we are doing in this work. 
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(a) (b) 



Figura 17: (a) Object space and (b) Image space for k = —6. 



We start with the analysis of the change of the image topology: in figurefTWa) we 
show the object (line segment given by Eq. (I48p with k = —6) and the caustic, which 
are in the plane z = zq = 11000. In figurefTWb) we show its image, which is not a 
segment of line but it is a distortion of it. For this value of k, that is for k = —6, we 
obtain only one image because any point of the line segment can be reached by only 
one reflected light ray. The first change in the image topology is obtained when the 
line segment is tangent to that part of the caustic corresponding to the circle of radius 
Rc, that is when k = —Rc, at this value of A; a new image is born as an isolated point. 
The second image that has born is visible for k = —4, as can be seen in figuresfTSTa) 
and llSlfb). In figure (jl9p we have presented the case k = —2. The following change 
of topology occurs when A; = 0, it is because the center of the circle with radius Rc 
is also part of the caustic. In figures [20r a)- l20r h) we show the object space and the 
image space for k = 0,2, 4, 6, respectively. Now it is clear that for this example, if 
k € {—Rc, 0) the object and its image have the same topology. The second image that 
was born when k = —Rc transforms into a circle when k = 0. For k € (0, Rc) we have 
two images. Finally, the last change in the image topology is obtained when k = Rc- 
For this value of k one of the images reduces to an isolated point, so that for k > we 
only have a single image. 

Remember that when k = —Rc in the image space there appears a new image from 
an isolated point, which as k goes to zero transforms into a circle, something similar 
happens for < A; < R^ but in this second case the circle reduces to an isolated point 
when k = Rc- In figure (I2ip we have isolated some of these loop images for some values 
of k. It is important to remark that this pattern was obtained by Ronchi [9] when 
he was studying a lens with a remarkable aberration. He also obtained the pattern 
corresponding to A; = 0. 

Finally, in the set of figures ([22]) and (f23]) we show the object space and the image 
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Figura 18: (a) Object space and (b) Image space for k = —4. 



space associated with a set of objects emulating the Ronchi ruhng and their associated 
ronchigrams. 




Figura 19: (a) Object space and (b) Image space for k = —2. 
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Figura 22: Simulation of the Ronchi ruling and the caustic curve. 
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Figura 23: Image of the Ronchi rulings or ronchigrams. 
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6. Conclusions 



In this work we have obtained an exact set of equations to study the change of 
image topology of an arbitrary object obtained by reflection on an arbitrary smooth 
surface. We have pointed out that the equations that allow to compute the image of an 
arbitrary one-dimensional object have been reported in the literature [13] in the context 
of the Ronchi and Hartmann tests. However, to our knowledge, nobody has remarked 
that the closed loops fringes observed in the ronchigram when the Ronchi ruling is 
located at the caustic place are due a disruption of images or fringes. We believe that 
the main contribution of this work is to realize that the procedure developed by Berry 
to study the changes of image topology of the Sun disk, is equivalent to that used in 
the simulation of ideal ronchigrams. From this observation it is clear that the caustic 
plays a major role in describing correctly this kind of pattern. We claim that the results 
established in this work could provide a description of the patterns obtained in other 
tests that use gratings. 

Our general results were illustrated by several examples when the reflecting surface 
is a spherical mirror and the point light source (in the Ronchi test) or observing eye (in 
Berry's procedure), is located on the optical axis. We believe, it could be worthwhile 
to study the change of image topology when the point light source is out of the optical 
axis and the reflecting surface is a conic reflector. 

We remark that analogous results can be obtained for lenses. In a future paper we 
will report these results. 
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